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Abstract: [ Aim] The aim of this study is to investigate the role of DaFOXO1 in the regulation of superoxide 
dismutase (SOD) gene expression and the developmental duration of summer diapause pupae of the onion fly, 
Delia antiqua. [ Methods] Genes encoding copper-zinc superoxide dismutase (DaCu/Zn SOD) and manganese 
superoxide dismutase (DaMn SOD) , downstream genes of DaFOXO1, were identified from the transcriptome 
data of D. antiqua. The features of amino acid sequences, subcellular localization and phylogeny of DaCu/Zn 
SOD and DaMn SOD of D. antiqua were analyzed by using bioinformatic tools. The expression patterns of 
DaFOXO1, DaCu/Zn SOD, and DaMn SOD at different developmental stages of summer diapause pupae of D. 
antiqua were determined by qRT-PCR. The effect of knocking down DaFOXO1 by RNAi on the expressions of 
DaCu/Zn SOD and DaMn SOD, the activity variations of DaCu/Zn SOD and DaMn SOD and the developmental 
duration of summer diapause pupae of D. antiqua were further analyzed. [ Results] The open reading frame 
(ORF) of DaCu/Zn SOD (GenBank accession no. KRO72551) of D. antiqua is 459 bp in length and encodes 
153 amino acids with a predicted molecular weight (MW) of 22.4 kD and an ioselectric point (pl) of 6.44, 
which belongs to a cytoplasmic Cu/Zn SOD. The ORF of DaMn SOD (GenBank accession no. KRO72549) is 
648 bp in length and encodes 216 amino acids with a predicted MW of 24.4 kD and a pl of 8.85, which 
belongs to mitochondrial Mn SOD. Amino acid sequence alignment revealed that DaCu/Zn SOD and DaMn SOD 
share 75% -94% identity with their homologues from other 10 species of Diptera, and contain typical SOD 
family domains. Phylogenetic analysis showed that they formed a robust phylogenetic branch with their 
homologues in Lucilia cuprina. qRT-PCR analysis revealed that the expression levels of DaFOXO1 were higher 
at the pre-diapause and post-diapause stages, but lower at the diapause stage. Higher expression of DaCu/Zn 
occurred at the diapause stage and post-diapause stage. However, the highest expression level of DaMn SOD 
was detected at the pre-diapause stage and diapause stage, and followed by at post-diapause stage. Knockdown 
of DaFOXO1 by RNAi significantly decreased the expression levels of DaCu/Zn SOD and DaMn SOD, and the 
activities of their corresponding enzymes, leading to a significantly extended duration of summer diapause 
pupae. [Conclusion] The results suggest that DaCu/Zn SOD and DaMn SOD are important members of FOXO1 
signaling network, and DaFOXO1 plays an important role in the regulation of the duration of summer diapause 
pupae of D. antiqua. 
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1 INTRODUCTION 


The onion fly, Delia antiqua, is an economical 
important pest of Allium plants throughout the 
northern hemisphere and can enter into summer or 
winter diapause to cope with unfavorable seasons 
(Hao et al., 2016). Insect diapause is a precisely 
regulated state that consists of several successive 
phases: induction, preparation, initiation, maintenance , 
termination and post-diapause quiescence ( Kostdl, 
2006). As an important adaption strategy, insect 
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diapause is thought to play important roles in 
preserving population and maintaining the population 
growth (Hao et al., 2016), which is characterized 
by arrested growth or development, decreased 
metabolism, increased stress resistance and extended 
lifespan ( Xu et al., 2018). Once diapause 
terminates, the development resumes at the post- 
diapause stage if conditions are favorable. However, 
clear molecular mechanisms of diapause termination 
are unknown, down-stream 


mechanisms involved in diapause hormone and 


especially the 
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ecdysone，and up-stream factors of termination. 
Therefore ，deeper understanding of the molecular 
mechanisms governing diapause termination is of 
of pest 
answers to 


high importance in the improvement 
management and also can give 
fundamental questions such as aging and lifespan. 
Previous studies revealed that many genes 
responsible for lifespan extension are conserved 


insulin/fork head 
transcription factor (INS/FOXO ) signaling is a key 


across phyla. Among them, 
regulator controlling the lifespan of diapause fly 
Drosophila melanogaster (Giannakou and Partridge, 
2007) , mosquito Culex pipiens (Sim and Denlinger, 
2008; Sim et al., 2015), 
Helicoverpa armigera (Zhang et al., 2017) , and the 


cotton bollworm 


dauer nematode Caenorhabditis elegans (Lee et al., 
2001 ). 
INS/FOXO signaling pathway is suppressed and 


During the whole diapause process, the 


regulates the expression of some genes associated 
with cold tolerance, UV resistance and reactive 
oxygen species (ROS) resistance (McElwee et al., 
2003; Murphy et al., 2003; Oh et al., 2006; 
Gershman et al., 2007). 

Several theories were put forward to address the 
underlying mechanisms of aging, including ROS 
detoxification mechanism, which is one of the well- 
known mechanisms for maintaining the homeostasis 


and lifespan extension. ROS includes O,, H,0O,, 
and OH’ , which can be generated as a byproduct of 
the normal metabolism. Excessive ROS can damage 
proteins, lipids and nucleic acids, and cause 
mitochondrial dysfunction, as well as serious effects 
on the integrity of cell membranes. To balance ROS 
levels, insects possess an effective antioxidant 


defense system that consists of enzymatic 


components, such as superoxide dismutase (SOD) , 
catalase (CAT) and peroxidase (POX), and non- 
enzymatic components, such as reduced glutathione 
and vitamin C (Jena et al., 2013). Among these 
antioxidant enzymes, SOD is the first and the most 
important enzyme against ROS. This protein family 
contains three major isoforms; cytoplasmic Cu/Zn 
SOD, extracellular Cu/Zn SOD, and mitochondrial 
Mn SOD. They differ in the metal binding, cellular 
location, structure and primary functions. However, 
the physiological level of ROS may act as messenger 
molecules in a variety of biological processes 
including health-extending ( Brys et al., 2010) and 
health-promoting ( Ristow et al., 2009). In dauer 
state of C. elegans, an equivalent stage of insect 
diapause, ROS can extend the 
inactivating the insulin signaling pathway ( Zarse et 
al., 2012). A recent study reported that the lifespan 


lifespan by 


of H. armigera diapause pupae was increased by the 
physical level of ROS (Zhang et al., 2017). These 
results indicate that insulin and ROS play important 
roles in diapause regulation. However, it is still not 
well-known how ROS plays a dual function and what 
signaling pathway is activated by ROS to extend 
lifespan of diapause insects. 

In our previous study, the expression levels of 
fork head transcription factor 1 gene ( DaFOXO1 ) 
and its putative downstream target gene DaCu/Zn 
SOD and DaMn SOD were differentially expressed 
during the diapause process of D. antiqua (Hao et 
al., 2016; Ren et al., 2018; Xu et al., 2018). 
Furthermore, the termination of diapause pupae of 
D. antiqua has not yet been investigated at the 
molecular level. Therefore, we asked whether 
FOXO-SOD signal is involved in the regulation of 
developmental duration of summer diapause pupae. 
To this end, downstream genes encoding DaCu/Zn 
SOD and DaMn SOD of DaFOXO1 of D. antiqua 
were identified from the transcriptome, and 
characterized by bioinformatic analysis. Moreover, 
the effect of knocking down DaFOXO1 by RNAi on 
the expression of DaCu/Zn SOD and DaMn SOD, 
their enzymatic activity 
duration were analyzed. Our results will be helpful 


for better INS/FOXO 
signaling pathway regulates the pupal duration of 


changes and diapause 


understanding how the 


onion fly and other diapause insects. 


2 MATERIALS AND METHODS 


2.1 Test insects 

A colony of D. antiqua was established in our 
laboratory ( Xu et al., 2018). Adults were 
maintained on artificial diet at 20 +0.5°C and 50% —- 
70% relative humidity with a photoperiod of 16L: 
8D. To obtain summer diapause pupae, adults and 
larvae were kept at 25 + 0.5°C and 50% - 70% 
relative humidity with a photoperiod of 16L: 8D. The 
pupation day was set as day 0 (DO). Pupae in pre- 
diapause (D0. 5 and D1. 5), diapause (D2, D4 
and D8) and post-diapause stages ( D14 and D20) 
were collected for survey of developmental stage- 
specific gene expression. Based on gene expression 
pattern, pupae at late diapause stage (D10) were 
selected for subsequent analysis. 
2.2 Identification of superoxide dismutase 
genes and bioinformatics analysis 

To identify genes 
BLASTP 


against the transcriptome data of onion fly using D. 


encoding superoxide 


dismutases , searching was performed 


melanogaster Mn SOD ( GenBank accession no. 


NP476925) and Cu/Zn SOD ( GenBank accession 
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no. NP476735) as queries ( Hao et al., 2016). The 
complete open reading frames ( ORFs ) were 
predicted using the ORF Finder (https: // www. 
ncbi. nlm. nih. gov/orffinder/). Conserved domains 
were predicted by the SMART (http; //smart. embl- 
heidelberg. de/) and CD search (https: // www. 
nebi. nlm. nih. gov/Structure/cdd/ wrpsb. cgi). The 
protein molecular weight (MW) and the isoelectric 
point (pl) were predicted by Expasy (https; // web. 
expasy. org/compute _ pi/). Signal peptide and 
subcellular localization were analyzed by SignalP3. 0 
(http; // www. cbs. dtu. dk/services/SignalP/) and 
WolfPsort ( https: // wolfpsort. hge. jp/), 
respectively. To know the phylogenic relationship of 
DaMn SOD or DaCu/Zn SOD with their homologues 
in closely related insect species, amino acid 
sequences were aligned using MAFFT (http: // 
mafft. cbre. jp/alignment/server/) , and ambiguous 
regions were removed by Gblock (http; // molevol. 
cmima. csic. es/castresana/Gblocks _ server. html ) . 
The best-fit evolutionary model was predicted by 
Modeltest ( Posada, 2006 ). The 
likelihood phylogenetic trees were constructed based 
on the model of WAG +G for both DaMn SOD and 
DaCu/Zn SOD with MEGA X ( Kumar et al., 


2018). The bootstrap consensus trees were inferred 


maximum- 


from 1 000 replicates. 
2.3 RNA extraction and cDNA synthesis 

Total RNA was extracted from one pupa of D. 
antiqua from Section 2. 1 using the TRIzol Reagent 
( Invitrogen, USA ) and treated with DNase to 
remove contaminated DNA. The quantity and quality 
NanoDrop ND-1000 


( Thermo Fisher Scientific) and gel electrophoresis, 


were determined with a 


cDNA was 
synthesized in a 20 uL reaction mixture using the 
RevertAid™ First Strand cDNA Synthesis Kit 
(Thermo, USA), including 3 wg RNA, 1 uL oligo 
d(T) 18 primer (50 ng/L) and incubated at 65°C 
for 5 min. Then, 5 uL 5 x reaction buffer, 1 uL 
Riolock™ RNase Inhibitor (20 U/L), 2 pL dNTP 
Mix (10 mmol/each), 1 pL RevertAid™ M-MuLV 
transcriptase (20 U/L) were added. The total 
reaction mix was incubated at 25°C for 5 min, 42°C 
for 60 min, 70°C for 5 min and stored at — 20°C for 
further use. 
2.4 qRT-PCR detection of gene expression 
patterns 

qRT-PCR was performed to detect the 
expression levels of DaFOXO1, DaCu/Zn SOD and 
DaMn SOD in D. antiqua at different developmental 
stages described as in Section 2. 1 in a 20 pL 


respectively. The first-stranded 


system, including 1 ML cDNA ( obtained from 
Section 2.3), 10 uL 2 x SybrGreen Mix, 0.8 pL 
each forward and reverse primer (10 ymol/L each) 
and 7.4 uL ddH,O. The PCR conditions were as 
follows: 94°C for 3 min; then 40 cycles of 94°C. for 
20 s, 58. for 20 s, 72°C for 20 s; and a final 
extension at 72°C for 3 min. The relative expression 
levels were calculated by the 2~“*" method (Livak 
and Schmittgen, 2001 ) using GAPDH as the 
reference gene. All samples were analyzed in three 
biological and three technical replicates. Gene 


specific primers were designed based on_ the 
sequences obtained in this study and the sequence of 
DaFOXO1 deposited in GenBank database 
( GenBank accession no. MG813258). All primer 


sequences were listed in Table 1. 


Table 1 Primers used in this study 








Gene GenBank accession no. Forward primer sequence (5' —3') Reverse primer sequence (和 —3’) Purpose 

DaFOXO1 MG8 13258 ACGCTTAATGGTCCGCTGTC TITGGATAGGCGGGTGACAT 

EGFP U55762 GCCACAAGTTCAGCGTGTC CTCGATGTTGTGGCGGATCT 

DaCu/Zn SOD KRO72551 CGTCATGTTGGTGATTTGGGT TCTTCGAATCGGGATGGGTT qRT-PCR 

DaMn SOD KRO72549 GTTGGGGTTGGTTGGGCTAT TCAACATAAGAGGGGCGCAA 

GAPDH NP525108 ACGTGGTGCTGCCCAAAACATCATT GGCGGACAGTCAAATCAACAAGG 

DaFOXO1 MG8 13258 CGTTCATTAACCTGCACCT GCCTGCAATTGATITATTTG 
TAATACGACTCACTATAGGC TAATACGACTCACTATAGGC 

人 MGS13258 GCGTTCATTAACCTGCACCT GGCCTGCAATTGATTTATTTG RNAi 

EGFP U55762 TAATACGACTCACTATAGGC TAATACGACTCACTATAGGCG 


GGCCACAAGTTCAGCGTGTC 


2.5 Preparation of double-stranded RNA 

A 596 bp fragment of DaFOXO1 was amplified 
from a previously constructed vector pMD-19T- 
DaFOXO1 using the primer set DaFOXO1-2 (Table 
1) as described by Xu et al. (2018). Template 
preparation for double-stranded RNA ( dsRNA ) 
synthesis was conducted in a 15 pL system 
containing 1 pL DNA template (596 bp fragment) , 


CTCGATGTTGTGGCGGATCT 


1.5 uL reaction buffer, 1.2 pL Mg (2.5 mmol/ 
L), 1.2 uL dNTPs (10 mmol/L each), 0.6 uL 
forward and reverse primer (10 mmol/L each ) 
(DaFOXO1-3), 0. 1 pL Taq polymerase and 8. 8 
uL H,O. The PCR program was as follows; 
denaturation at 94°C for 3 min; 35 cycles of 94°C for 
30 s, 56 for 30 s, and 72°C for 1 min; and a final 
extension at 72°C for 3 min. The PCR product was 
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checked in a 1.0% agarose gel, and then purified 
for dsRNA preparation. The dsRNA synthesis was 
conducted in a 20 uL reaction mixture, including 10 
uL RiboMAX™ Express T7 2 x buffer, 1 uL above 
PCR product, 7 uL H,O and 2 uL enzyme mix. The 
mixture was incubated at 37°C for 12 h, and then 
1 uL RNase-free DNase (1 U/uL) was added to 
remove DNA, followed by the incubation at 95°C for 
5 min, 72°C for 10 min, 60°C for 10 min, 50°C for 
8 min and 37° for 10 min. The dsEGFP was also 
prepared and used as the positive control as 
described previously ( Xu et al., 2018). The 
resulting synthesized dsDaFOXO1 or dsEGFP was 
precipitated with LiCl solution, suspended in RNase- 
free water and quantified by a NanoDrop ND-1000 
(Thermo Fisher Scientific). The concentration was 
quantified and diluted to a final concentration of 
1 mg/ WL for later use. 
2.6 Effect of knocking down DaFOXO1 by 
RNAi on the expression levels of DaMn SOD and 
DaCu/Zn SOD 

Three hundred and sixty summer diapause 
( D10 ) pupae were divided into three groups 
randomly (120 pupae each). Based on the dose 
optimization, each pupa was injected with 1 pg 
dsDaFOXO1 (treatment group) or dsEGFP ( positive 
After the 
injection, all pupae were transferred to the original 


control ) via the abdomen region. 
rearing condition (see Section 2. 1). The pupae 
without any treatment were used as the negative 
control (CK). The sampling time was set at 0, 4, 
8, 12, 16, 20, 24 and 36 h after injection and three 
pupae at each sampling point for each group were 
collected for gene expression pattern analysis by 
qRT-PCR as described in Section 2. 4. 
2.7 Effect of knocking down DaFOXO1 by 
RNAi on SOD activities 

RNAi treated pupa from Section 2.6 (3 pupae 
for each sampling point) was homogenized in 0. 5 
mL physical salt solution, centrifuged at 4°C 3 000 
r/min for 5 min and the supernatant was used for 
protein quantification and enzymatic activity assays. 
The total protein content was determined using a 
commercial kit ( Comin Biotech, Suzhou). The 
activities of DaMn SOD and DaCu/Zn SOD were 
determined using the Total SOD (T-SOD) detecting 
Kit ( Nanjing Jiancheng Bio. Nanjing ) 
following the manufacturer’ s instructions. The 
DaCu/Zn SOD activity was obtained by subtracting 
the activity of DaMn SOD with inhibitor addition. 
The SOD activity was expressed in unit ( U) per 


Inst., 


microgram total protein. 


2.8 Effect of knocking down DaFOXO1 by 
RNAi on the pupal eclosion 

To investigate the effect of knocking down 
DaFOXO1 by RNAi on the pupal eclosion, 90 pupae 
for each group from Section 2. 6 were selected for the 
eclosion experiment. Newly eclosed flies were 
recorded daily and the eclosion curve for adult 
emergence was plotted with Excel 2010. 
2.9 Statistical analysis 

All data were presented as means + SD of three 
biological and three technical replicates except the 
pupal eclosion experiment. One-way analysis of 
variance ( ANOVA ) was performed for time-course 
expression pattern analysis and followed by a 
Duncan’ s test. Significant differences among data 
from three groups at the same sampling time were 
determined using a Student’ s t-test. All statistical 
analyses were performed using SPSS19. 0 ( SPSS, 
Inc., Chicago, IL, USA). The statistical difference 
was considered as significant at P<0.05 and 
extremely significant at P <0. 01. 


3 RESULTS 


3.1 Gene identification, characterization and 
phylogenetic analysis 

In this study, DaMn SOD ( GenBank accession 
no. KRO72549 ) and DaCu/Zn SOD ( GenBank 
accession no. KR072551 ) were identified based on 
the transcriptome data. The ORF of DaCu/Zn SOD 
is 459 bp in length, coding a protein of 153 amino 
acids with a predicted MW of 22.4 kD and a pl of 
6.44. Subcellular localization prediction revealed 
that DaCu/Zn SOD is a cytoplasmic protein. 
Similarity analysis showed that DaCu/Zn SOD shares 
94% amino acid sequence identity with LeCu/Zn 
SOD (GenBank accession no. XP023303945 ) of 
Lucilla cuprina, 91% identity with MdCu/Zn SOD 
(GenBank accession no. NP001295981) of Musca 
domestica and 88% identity with ScCu/Zn SOD 
( GenBank XP013106323 ) of 
Stomoxys calcitrans. Multiple sequence alignment 
revealed that DaCu/Zn SOD contains the conserved 
domain of Cu/Zn SOD family (a. a. 22 - 176, E = 
le-49) (Fig. 1; A). Cu’*-binding sites (H*, 
H”, H? and H'”) and Zn’* -binding sites ( H®, 
H”, H”, and D”) are also conserved (Fig. 1; A). 
The phylogenetic tree showed that DaCu/Zn SOD 
was phylogentically clustered with LcCu/Zn SOD of 
L. cuprina (Fig. 1; B). 

The ORF of DaMn SOD is 648 bp in length, 
encoding a 216-amino-acid protein with a predicted 
molecular mass of 24. 4 kD and a pl of 8. 85. 
Subcellular localization results showed that DaMn 


accession no. 
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Lucilia cuprina LeCu/Zn SOD (XP 023303945) 
Delia antiqua DaCu/Zn SOD (ALF39391) 
Musca domestica MdCu/Zn SOD (NP001295981) 

Stomoxys calcitrans ScCu/Zn SOD (XP013 106323) 
Rhagoletis zephyria RzCu/Zn SOD (XP017470220) 


Bactrocera dorsalis BdCu/Zn SOD (NP001291905) 


Zeugodacus cucurbitae ZcCu/Zn SOD (XP011176600) 


Ceratitis capitata CeCu/Zn SOD (XP004522289) 


Drosophila melanogaster DmCu/Zn SOD (NP476735) 


100 Drosophila ananassae DraCu/Zn SOD (XP001956494 ) 


m Drosophila bipectinata DbCu/Zn SOD (XP017090148) 


— 
0.050 


Fig. 1 


Multiple sequence alignment (A) and phylogenetic tree (B) of DaCu/Zn SOD of Delia antiqua 


with its homologues from other 10 species of Diptera by maximum likelihood method 


In Fig. A, the consensus head and tail were boxed, the active sites responsible for Cu2+ - or Zn2+ -binding sites were marked with inverted triangles, and 


the signatures of Cu/Zn SOD family were underlined. In Fig. B, the phylogenetic tree was constructed under the best evolutionary model WAG + G. 


Bootstrap values from 1 000 replicates were shown on each branch. The scale bar represents the substitution per amino acid site. 


SOD is localized in the mitochondria. It shares 90% 
amino acid sequence identity with mitochondrial MdMn 
SOD ( GenBank accession no. NP001273816) of M. 
dometica and 92% identity with mitochondrial LeMn 
SOD (GenBank accession no. XP023307201 ) of L. 
cuprina. DaMn SOD contains an N-terminal domain 
(a.a.19 —- 100, E =7.3e —34) and a C-terminal 
domain (a. a 104 - 211, E =9e - 37) (Fig. 2: 
A). In addition, DaMn SOD contains a typical Mn 
SOD signature ( DVWEHAYY ) and four Mn’*- 
binding sites ( H“, H®, D'” and H"™'). The 
phylogenetic analysis revealed that DaMn SOD was 
grouped with Mn SOD of L. cuprina, M. domestica 
and S. calcitrans (Fig. 2: B). 
3. 2 Developmental expression analysis of 
DaFOXO1, DaCu/Zn SOD and DaMn SOD in 
summer diapause pupae 

To gain insights into the importance of FOXO- 
SOD signal in summer diapause pupae, gene 


expression patterns in different developmental stages 
of D. antiqua were determined by qRT-PCR. The 
results showed that the transcription level of 
DaFOXO!1 was suppressed at the diapause stage (D2 — 
D8) when compared with those at the pre-diapause 
stage (DO.5 - D1. 5), but the highest expression 
level occurred at the post-diapause stage ( D20 ) 
(Fig. 3: A). Higher expression levels of DaCu/Zn 
SOD was detected at the diapause stage (D2) and 
post-diapause stage (D20) (Fig. 3: B). However, 
no clear expression pattern was found for DaCu/Zn 
SOD. The highest expression level of DaMn SOD 
was found at the pre-diapause stage (D1.5) and 
diapause stage (D4), then followed by the post- 
diapause stage ( D20) (Fig. 3: C). 

3.3 Effect of knocking down DaFOXO1 by 
RNAi on the expression of DaCu/Zn SOD and 
DaMn SOD 


The results showed that the transcription levels 
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Bactrocera dorsalis BdMn SOD (XP01 1203981) 
Bactrocera latifrons BIMn SOD (XP018788990) 


Zeugodacus cucurbitae ZcMn SOD (XP011184151) 


Rhagoletis zephyria RzMn SOD (XP017478302) 
Ceratitis capitata CeMn SOD (XP004517465) 
Drosophila guanche DgMn SOD (SPP72797) 


Drosophila persimilis DpMn SOD (XP026841012) 


Delia antiqua DaMn SOD (ALF39393) 


Lucilia cuprina LeMn SOD (XP023307201) 
Musca domestica MdMn SOD (NP001273816) 


Stomoxys calcitrans ScMn SOD (XP013109790) 


Fig. 2 Multiple sequence alignment (A) and phylogenetic tree (B) of DaMn SOD of Delia antiqua and 
its homologues from other 10 species of Diptera by maximum likelihood method 
In Fig. A, the mitochondrial transit peptide, the Tom20 recognition motif and the signature of Mn SOD family were underlined, and the active sites 


involved in Mn?* 


- binding were marked with inverted triangles. In Fig. B, 


the phylogenetic tree was constructed under the best evolutionary model 


WAG +G. Bootstrap values from 1 000 replicates were shown on each branch. The scale bar represents the substitution per amino acid site. 


of DaCu/Zn SOD and DaMn SOD were not affected 
by the injection of dsEGFP (Fig. 4: B, C). 
However, the expression levels of DaCu/Zn SOD 
(Fig. 4: B) and DaMn SOD (Fig. 4: C) were 
significantly decreased after injection of dsDaFOXO1 
(all P <0. 05) as compared with the negative 
control (CK) and the suppression effect lasted for 36 
h. The minimum expression levels of DaCu/Zn SOD 
and DaMn SOD were observed at 8 and 4 h after 
injection of dsDaFOX01, respectively (Fig. 4). 
3.4 Effect of knocking down DaFOXO1 by 
RNAi on the activities of DaCu/Zn SOD and 
DaMn SOD 

To investigate whether knocking down 
DaFOXO1 contributed to the reduction of SOD 
activities, variations were 
measured at 0 - 36 h after the microinjection. The 


enzymatic activity 


results showed that both activities of DaCu/Zn SOD 
(Fig. 5: A) and DaMn SOD (Fig. 5: B) were 
suppressed from 4 h after injection of dsDaFOXO1 
when compared to the dsEGFP-injected group and 
the negative control (CK) at each sampling time 
point (all P<0.05). The minimum activities of 
DaCu/Zn SOD and DaMn SOD were observed at 24 
and 8 h after injection, respectively. There was a 
clear time lag between variations of gene expression 
and the corresponding enzyme activity. 
3.5 Effect of knocking down DaFOXO1 by 
RNAi on the pupal duration 

To elucidate the function of FOXO-SOD signal 
in the maintenance phase of summer diapause, the 
effect of knocking down DaFOXO1 on the pupal 
duration was determined based on the eclosion curve 


of diapause pupae. Interestingly, the length of 
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Fig. 3 Expression patterns of DaFOXO1 (A), DaCu/Zn SOD (B) and DaMn SOD (C) 


in summer diapause pupae of Delia antiqua 
Data in the figures are mean + SD of three biological replicates and three technical replicates for each sampling time point (n =3 x3). Gene expression 
levels in each group were evaluated using one-way analysis of variance (ANOVA) followed by a Duncan’ s multiple comparison test, and different letters 
above bars indicate significant difference ( P <0.05). 
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Fig. 4 Changes in relative expression levels of DaFOXO1 (A), DaCu/Zn SOD (B) and DaMn SOD (C) in 
summer diapause pupae of Delia antiqua after RNAi of DaFOXO1 detected by qRT-PCR 
Pupae (D10) injected with dsEGFP were used as the positive control, while those without any treatment were taken as the negative control (CK). The 
house-keeping gene GAPDH was used as an internal reference gene for data normalization. The same for Figs. 5 - 6. The expression level of each gene in 
samples at 0 h (just after injection) was used as the calibrator. All data are mean + SD of three biological replicates and three technical replicates (n = 
3 x3). Pairwise comparisons of the difference within each sampling time point were performed by the independent samples t-test. The difference was 


considered as significant when * P <0.05 or “ P<0.01. n.s. represents no significant difference. The same for Fig. 5. 


diapause maintenance phase was significantly (CK). The eclosion peak was postponed about 6 
extended in both dsEGFP and dsDaFOXO1 injected and 15 d in dsEGFP and dsDaFOXO1 injected 


groups when compared with the negative control groups, respectively (Fig. 6). 
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Fig.5 Effect of knocking down DaFOXO1 by RNAi on the activities of DaCu/Zn SOD (A) and 
DaMn SOD (B) in summer diapause pupae of Delia antiqua 
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Fig. 6 Effect of knocking down DaFOXO1 by RNAi on the 
duration of summer diapause pupae of Delia antiqua 
Pupae ( D10) injected with dsDaFOXO1 or dsEGFP were kept at the 
same condition as the negative control (CK). Eclosion events were 
recorded daily over 40 d. Cumulative numbers of eclosed individuals 


were plotted against eclosion day using EXCEL 2010. 


4 DISCUSSION 


Some insects can enter into diapause, a 
phylogenetically regulated state of dormancy, which 
is characterized by the reduced metabolism, 
developmental arrest, and high stress tolerance. This 
terminated by favorable 
increasing or 


arrest state can be 
environmental cues ( such as 
decreasing temperature, photoperiod change ) or 
chemical treatment ( like hormone stimuli ) 
( Ragland et al., 2011 ). The  post-diapause 
development resumes if conditions are favorable, 
which is a very complicated process involved in 
various cellular events, such as protein synthesis, 
increased metabolism, cellular differentiation and 
morphological changes ( Nambu et al., 1997; 
Malarkey et al., 2010). All these events occur 
together with large energy fluctuations and promote 
the accumulation of reactive oxygen species (ROS) , 
such as superoxide, and 
hydroxyl radicals. Excessive ROS generated as by- 
products damage cellular 
components. However, physiological ROS level can 


hydrogen peroxide, 


of metabolism can 


participate in cellular signaling transduction and is 
essential for lifespan extension (Finkel, 2003). The 
involvement of ROS in cellular signaling is part of a 
physiological regulation. 


complex, redox-based, 


Under the normal condition, a delicate balance of 
intracellular oxidative stress produced by ROS is 
achieved by a complex system of antioxidant defense 
(Hao et al., 2018). Superoxide dismutase enzymes 
(SOD) are the most important regulators of ROS 
level ( Landis and Tower, 2005). Cu/Zn SOD 
functions in the cytoplasm or outer mitochondrial 
space, while Mn SOD functions in the inner 
mitochondrial space. In Drosophila, the mutation of 
Cu/Zn SOD or Mn SOD reduced lifespan ( Phillips 
et al., 1989; Kirby et al., 2002) , while their over- 
expression could extend the longevity ( Parkes et al., 
1998; Sun et al., 2004). 

Previous study showed that the suppression of 
insulin signaling pathway gene could extend the 
lifespan of dauer larvae of C. elegans and Drosophila 
adults. It was thought that this process was regulated 
by the insulin-target transcription factor FOXO 
(Clancy et al., 2001; Tatar et al., 2001). Two 
SOD genes were identified in C. pipiens: a 
mitochondrial Mn SOD gene ( sod-1 ) 
cytoplasmic Cu/Zn SOD gene (sod-2). However, 
only sod-2 was under the control of FOXO and 
protected ovaries of winter diapause adults (Sim and 
Denlinger, 2011). Why sod-2 is the only SOD gene 
involved in the diapause response of C. pipiens is not 


and a 


immediately clear. Interestingly, the expressions of 
orthologues of sod-2 and sod-3 were significantly up 
regulated and involved in lifespan regulation of C. 
elegans (Murphy et al., 2003). 

In our previous work, the involvement of 
antioxidant enzymes in redox level balance during 
the whole summer or winter diapause process had 
been determined ( Hao et al., 2018). Moreover, we 
have found that the insulin/FOXO signaling pathway 
plays an important role in the induction and 
maintaining of summer diapause in D. antiqua (Xu 
et al., 2018). However, the regulatory role of 
FOXO-SOD signal in lifespan or summer diapause 
maintenance of onion fly is not known. In this study, 
two SOD genes were identified based on the 
transcriptome data and their proteins were putatively 
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localized in cytoplasm ( DaCu/Zn SOD ) and 
mitochondria (DaMn SOD). Their expression levels 
were coordinated with the expression of DaFOXO1. 
Knockdown of DaFOXO1 by RNA interference can 
suppress the expressions of DaCu/Zn SOD and 
DaMn SOD at both transcriptional and translational 
levels, but the suppression degree was different. 
Their different temporal and spatial expression 
patterns suggest that they contribute to antioxidant 
protection at different developmental periods. Low 
of DaFOXO\ and DaSODs can 
significantly extend the diapause duration of summer 
diapause pupae. 
specific expression pattern of these genes may help 
explain how FOXO-SOD 


development of diapause 


expression 
Further examination of tissue- 
signal regulates the 
pupae. 
diapause pupae injected with buffer with dsEGFP 


longer diapause Although 
evidence indicated that water was often considered as 


Interestingly , 


also have duration. 
essential for diapause complete and prerequisite for 
the resumption of development ( Hodek , 1996) , it is 
argued that this phenomenon is caused by water or 
dsEGFP. Similarly, when the larvae of the Colorado 
beetle, Leptinotarsa decemlineata, were fed with 
dsEGFP, the period was 


extended, which causes the delay of the pupation 


larval significantly 


(Guo, 2015). But the mechanism was unknown. 
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E, [RK] AAR K EAA B38 Delia antiqua & iit FAK A DaFOXO1 对 超 氧 化 物 上 化 酶 
(SOD) 基 因 表 达 及 肾 发 育 历 期 的 调控 作用 。【 方 法 】 从 草 晶 转录 组 数据 中 鉴定 DaFOXO1 下 游 铜 锌 
超 氧化 物 歧化 酶 基因 DaCu/Zn SOD 和 锰 超 氧化 物 歧 化 酶 基因 DaMn SOD; 利 用 生物 信息 学 工具 对 
DaCu/Zn SOD 和 DaMn SOD 的 和 氨基酸 序 列 特征 、 亚 细胞 定位 和 系统 发 育 关 系 进行 分 析 。 通 过 qRT- 
PCR 方法 分 析 DaFOXO1, DaCu/Zn SOD 和 DaMn SOD 2B EYE AAR AR PY RIK 

特点 ;进一步 分 析 DaFOXO1 基因 被 干扰 后 , 中 DaCu/Zn SOD 和 DaMn SOD 基因 的 表 
AA ME EAS Be Le AAA Bm, (AR) BES) Bw DaCu/Zn SOD 
(GenBank 登录 号 : KRO72551 ) 的 开放 阅读 框 长 459 a 153 个 氨基 酸 ,预测 蛋白 分 子 量 为 
22.4 kD, 等 电 点 为 6.44, 属于 细胞 质 型 铜 锌 超 氧 化 此 化 酶 ; DaMn SOD ( GenBank 
KR072549 ) 的 开放 阅读 框 长 648 bp ,编码 216 个 和 氨基酸, 预测 蛋白 分 子 量 为 24.4 kD ,等 电 点 为 
8. 85 , 属于 线粒体 型 锰 超 氧化 物 歧 化 酶 。 和 氨基 酸 序列 比 对 结果 显示 ,DaCu/XZn SOD 和 DaMn SOD 与 

其 他 10 FRM ARMA RESA 5% ~94% 的 和 氨基 酸 序 列 一 致 性 , 且 具 有 典型 的 SOD 家 族 序 
列 特征 ; 系统 发 育 分 析 显 示 它 们 与 铜绿 蝇 Lucilia cuprina 同 源 蛋 和 白 形 成 高 支持 率 的 一 支 。qRT-PCR 
分 析 表 明 ,DaFOXO1 基因 在 滞 育 前 期 和 滞 育 后 期 的 表达 量 较 高 ,而 在 滞 育 期 的 表达 量 低 ; DaCu/ 
Zn SOD 基因 在 滞 育 期 和 滞 育 后 期 哇 高 表达 ;但 DaMn SOD 基因 在 滞 育 前 期 和 灌 育 期 的 表达 量 最 
高 ,在 滞 育 后 期 次 之 。 干扰 DaFOXO1 可 显著 抑制 DaCu/Zn SOD 和 DaMn SOD 的 基因 表达 及 相应 2 
酶 活性 ,并 能 明显 延长 夏 滞 育 肾 的 滞 育 期 。【 结论】 结果 说 明 , DaCu/Zn SOD 和 DaMn SOD 是 
FOXO1 信号 网 络 中 的 重要 成 员 ;DaFOXO1 sp B48 Ze A A A ERIE A 
关键 词 : Bx; FOXO1; WA: 超 气 化 物 歧化 酶 ; BA; RNAI 
中 图 分 类 号 : Q966 ”文献 标识 码 : A ”文章 编号 : 0454-6296(2019)12-1359-10 
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